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Abstract: Optical binding interactions between laser-trapped spherical 
microparticles are familiar in a wide range of trapping configurations. 
Recently it has been demonstrated that these experiments can be accurately 
modeled using Mie scattering or coupled dipole models. This can help 
confirm the physical phenomena underlying the inter-particle interactions, 
but does not necessarily develop a conceptual understanding of the effects 
that can lead to future predictions. Here we interpret results from a Mie scat- 
tering model to obtain a physical description which predict the behavior and 
trends for chains of trapped particles in Gaussian beam traps. In particular, 
it describes the non-uniform particle spacing and how it changes with the 
number of particles. We go further than simply demonstrating agreement, 
by showing that the mechanisms "hidden" within a mathematically and 
computationally demanding Mie scattering description can be explained in 
easily-understood terms. 
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1. Introduction 

The term optical binding describes the light-mediated lateral or longitudinal interactions be- 
tween groups of microparticles trapped by one or more laser beams QJ. Through their coherent 
scattering properties, particles will modify the light field surrounding them, and in this way 
multiple particles are able to interact through the medium of the electromagnetic field around 
them. This gives rise to a rich tapestry of nonlinear static and dynamic behavior, including chain 
formation J2j[3], bistability |4|, two-dimensional "crystal" arrays J5) and periodic particle mo- 
tion Ed. 

Attempts have been made to model such experiments using Mie scattering models J6] 111 HI 
and coupled dipole models J9). Extremely good agreement was recently reported between a 
Mie scattering model and the behavior of particles trapped in a counter-propagating Gaussian 
fiber-based trap |7|, and a model based on coupled dipole calculations has also recently been 
to inform a simple explanation of the mechanisms which lead to chain formation in a counter- 
propagating Bessel beam trap (9J . 

In this paper we will use numerical results from a Mie scattering model, as well as from a 
simpler heuristic model, to build up a detailed understanding of the mechanisms which lead 
to optical binding of a one-dimensional chain of trapped particles in a counter-propagating 
Gaussian beam trap. This was possibly the first configuration in which ID optically bound 
chains were observed Q. The inter-particle spacing decreases as the number of particles in the 
trap increases, and a slight anisotropy is observed in the chain whereby the inter-particle spacing 
towards the center of the chain is smaller than the spacing at the edges of the chain [7. Fig. 5]. 
Any comprehensive model of the optical interaction must be able to explain these results. 

We will show that despite the apparent similarities between the trapped chains in this config- 
uration and the trapped chains in counter-propagating Bessel beams J5]|, the binding mechanism 
is entirely different for a Gaussian beam trap, and is dominated by "radiation pressure" effects 
(the scattering force iflOl ). Now that it is clear that Mie scattering models can accurately re- 
produce experimental results first reported nearly seven years ago, the challenge is to interpret 
those results in easily-understood terms, and to distill out the key mechanisms "hidden" within 
the complex model in order to develop a conceptual understanding of why optical binding and 
chain formation occurs. We will explain how the trapped particles modify the beam shape to 
support stable chains of particles with non-uniform inter-particle spacings, whose spacing de- 
creases as the number of particles in the trap increases. 

2. Optical binding concepts, and modeling 

When light is incident on a particle, the particle scatters the light, producing a secondary field 
which radiates in all directions. An important question to be asked is: what governs the spacing 
of the particles - the forward-scattered field or the back-scattered field? As has been observed 
by numerous authors ifTTl [121 [131 [9), the interference of back-scattered light with the inci- 
dent forward-propagating beam leads to a very large number of nearby trapped configurations, 



separated by either one or half a wavelength depending on the geometry. In the Gaussian beam 
trapping geometry we are considering here, the typical particle spacings are considerably larger 
than the wavelength of the laser Q (and indeed the particles may be several wavelengths in di- 
ameter). Consequently this binding mechanism, which is most significant for particles up to 
about half a wavelength in diameter, has little effect on the large-scale behavior of the trapped 
chains, which is determined largely by the forward-scattered light. 

Because of this important distinction between the effects of forward- and back-scattering, all 
the models discussed in this paper will explicitly treat the two incoherent counter-propagating 
beams separately: it is important that effects due to forward-scattering of one beam can be sepa- 
rated from back-scattering effects due to the other beam if the mechanisms are to be understood. 
Within a theoretical model, we can not only separate the effects of forward- and back-scattering, 
but can even "disable" one of these effects. The Mie scattering model for N particles can be 
expressed in the following linear algebra form lfl4ll : 

s (,) = T.a« (1) 

where a' ! ' represents the net field incident on particle i, a^ t the zero-order incident field on 
particle i due to the external laser field in the absence of other particles and s^ the scattered 
field from particle i; F ; -, is the translation matrix from a basis centered on particle j to one 
centered on particle i, and T is the T-matrix representing the scattering properties of the particle, 
which is diagonal for spherical particles. This expression simply states that the field incident 
on a particle i is the (coherent) sum of the laser field and the field scattered by all the other 
particles. 

Using this framework, it is easy to alter the model so that backscatter is not taken into ac- 
count. If the particles are indexed "upstream" to "downstream" (i.e. with the first particle closest 
to the laser source), then we simply modify the sum to read: 

a«= a «+£F;,. S W 

Although such a model is un-physical, there can be significant benefits in modeling such a situa- 
tion: a real-world experiment will be subject to Brownian motion of the trapped particles, which 
will wash out small local minima in the interactions caused by the back-scattered field. In the 
absence of Brownian motion (i.e. at absolute zero), a simulation can easily become trapped in 
such a local minimum. Since it is computationally very expensive to include Brownian motion 
in a simulation, choosing a model which does not include the back-scattered light is a useful 
compromise in many cases (though the results of these un-physical simulations must be verified 
using a full physically -realistic model). Other details of the model used can be found in [H. 

3. Chain formation in Gaussian beam traps 

In order to understand the binding mechanisms, consider first the intensity distribution down- 
stream of a single 1.0 fim diameter particle illuminated by a single plane wave. The dominant 
effect is a "focusing" of the light (in the limit of large particles we can consider the particle as 
a spherical lens within the framework of ray optics). Consider now the force on a second par- 
ticle placed in this field. This is plotted in Figure Q] which shows how two particles are stably 
bound in a Gaussian beam trap, but are not stably bound in counter-propagating plane waves. 
We can see that the force due to the light focused by the first particle causes the two particles 
to be repelled in the case of counter-propagating plane waves (Figure [TJ)). When we consider 
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Fig. 1. Force on a particle in a two particle system, as a function of particle spacing, (a) 
Force on the downstream (right-hand) one of a pair of 1 flm particles illuminated with 
a single right-going plane wave; (b) repulsive force between the pair of particles when 
two counter-propagating beams are used; (c) as (a) but for a 25mW Gaussian beam; (d) 
two counter-propagating Gaussian beams, showing the stable spacing of ~ 8 flm. Note the 
modulation due to backscattered light, which has little effect on the general trend of the 
binding behavior. The broad harmonic potential introduced by the use of Gaussian beams 
has altered curve (b) to give curve (d), which is outwardly similar but which has a stable 
inter-particle spacing at around 7 flm (marked with an arrow). 



counter-propagating Gaussian beams, through symmetry there is no net force on the center of 
mass of the particle pair, and the beams provide a broad background harmonic trapping poten- 
tial. The particles will stabilize with a spacing which is largely determined by the balance of 
the repulsive force between the two particles and the harmonic trapping potential of the trap 
(Figure [TJl), as suggested in JSJ. We emphasize that although we refer to a "focusing" of the 
light, we are far from the ray-optics regime, and it is not appropriate to use a ray-optics formula 
for the focal length, or to suggest that one particle will be bound "at the focus" formed by the 
other particle. 

If we just consider the effects of a single beam, then in this experimental setup the contribu- 
tion of the gradient force turns out to be only a small fraction of the total force on a particle. 
However, remember that with two counter-propagating beams a large part of the force exerted 
by one beam is balanced by the force due to the other beam. With some experimental param- 
eters the contribution of the gradient force to the net force on a particle due to the two beams 
together can be non-negligible. Thus it is not really possible to state that either the gradient 
force or the scattering force will dominate under all circumstances. 

A more interesting case than the two-particle case is to consider is that of a larger number 
N of trapped particles (indexed i = 1 to N), for which we intend to explain the three main 
properties of the particle chains: 

• The fall in inter-particle spacing with N. 

• The anisotropic particle spacing within the chain, with a smaller particle spacing near the 
middle of the chain. 



• For some experimental parameters, chains are only supported up to a certain number of 
particles, beyond which the chains collapse. 

There are a number of statements we can make about the behavior of the chains, based on 
symmetry considerations, with little or no assumptions on the nature of the inter-particle inter- 
actions: 

1 . Since the arrangement of the beams (two incoherent beams counter-propagating along 
the z axis) is symmetric about the z = plane, the force on particle ; due to one beam is 
equal and opposite to the force on particle N—i+1 due to the other counter-propagating 
beam. This is true for any symmetric arrangement of particles, whether or not this is 
an equilibrium configuration. If / ; + (or fr~) is the force on particle i due to the beam 
propagating in the +z (or — z) direction, then = f^_ i+l - 

2. In addition, in equilibrium, there must be a net force of zero on each particle when the 
forces from the two beams are added together (i.e. / ; + = ff), since by definition there 
must be no particle motion in equilibrium. Combining this with the previous requirement, 
we have = f^_ i+l - In other words, the forces on the particles in the chain due to each 
individual beam must be symmetric about the center of the chain. 

In addition to the above, we will make a number of simplifying hypotheses about the interaction 
mechanism between the particles. These hypotheses have been observed to be approximately 
true in numerical experiments: 

3. The particle spacings are determined by/brwani-scattering. 

4. The force on a given particle i is a function of the light intensity /q which would be 
found be at that point in the absence of that given particle (the Born approximation; see 
figure |2}. 

5. The profile of the on-axis scattered intensity downstream of a given particle ; has the 
form /q ' X I(z — z/), where /(z — z/) is a fixed downstream intensity profile which applies 
to any particle at any position. Consequently, the force on particle i + 1 has the form 
/q x F(z — Zj) for a fixed downstream force profile F(z — Z/). 

These statements are already enough to explain why the inter-particle spacing at all points in 
the chain will decrease if an additional particle is added onto either end of the chain. The 
justification for this is as follows. The force pushing inwards on what is now particle 2 in the 
chain has been increased (it was previously just the force due to the unperturbed laser beam; 
it is now enhanced by the additional light focused onto it by particle 1). Assuming there are 
some losses along the length of the chain, then the force pushing outwards on what is now the 
last-but-one particle in the chain will also increase, but by a smaller amount. Hence whenever 
additional particles are added to the chain, the inner ones will be pushed inwards, and so the 
inter-particle spacing between any given pair of particles in the chain will decrease. Equilibrium 
is then restored because the closer inter-particle spacing enhances the transmission efficiency, 
thereby further increasing the repulsive force on the last particle in the chain. 

In order to explain the non-uniform particle spacing, we propose a simple ansatz model for 
the force fi as a function of z; and z,-i, as follows: 

fi(zuZi-i) = e-^ x /,_! + (/„ - j3z«) (2) 

Here the first term represents assumption|5]and the second term represents a background inten- 
sity due to the laser field, which is decreasing with distance from the beam waist. We emphasize 
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Fig. 2. A simple Born approximation model uses the intensity of a right-going Gaussian 
beam in the absence of particle 3 (whose location is indicated with an arrow in (c)) to 
determine the force f$ on that particle, (a) shows the field for 3 particles exposed to a 
single Gaussian beam; (b) shows the field in the absence of particle 3 and (c) plots this field 
(all data in this figure was generated using a full Mie scattering model). 



that the functional form of fi has simply been selected empirically to give a reasonable approx- 
imation to the observed inter-particle force. If a closer agreement with the Mie scattering model 
was desired, a "hybrid" model could be used, in which /, is actually determined from the inter- 
particle forces for a pair of particles in a plane wave, calculated using Mie scattering theory. 
However, we have instead chosen to keep our model as elementary as possible. 

Having made these assumptions, we can test the predictions of this very simple model against 
the definitive calculations of a rigorous Mie scattering model, and against established experi- 
mental observations. Figure[3]illustrates how this model predicts the force to vary along a chain 
of particles with a constant inter-particle spacing, and an example of how the forces calculated 
from a Mie scattering model vary along a similar chain with constant particle spacing. It can 
be seen from this plot that requirement |2] (a symmetric force profile) is not satisfied in either 
model: the force pushing particle 1 inwards is weaker than the force pushing particle N out- 
wards, which means that the tendency will be for the outermost particles to move apart. The 
natural next step is therefore to allow the particle spacings to vary along the length of the chain, 
as they would do in real life in response to this repulsive force. As Figure [4] shows, this ap- 
proach allows a symmetric force profile to be produced, and leading to stable trapping of the 
chain with these slightly non-equilibrium spacings. 

The first few particles in the chain act to focus the laser field onto the next particle in the 
chain, and hence initially the force rises sharply with particle index ;. Particles towards the 
middle of the chain can be thought of as acting more like a (very inefficient) waveguide where 
the intensity is propagated from one particle to the next with some losses, which are compen- 
sated for by the re-focusing of additional background light. The intensity (and force) then drops 
again towards the end of the chain due to the increased particle spacings. 
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Fig. 3. Example of how the force (arbitrary units) on particle i in a chain of 8 particles due 
only to the right-going beam varies along the chain if the inter-particle spacing is constant. 
Mie scattering model (red, unbroken line) and our simple ansatz model (green, dotted line). 
Both plots show similar trends; neither is symmetric with respect to the center of the chain. 
This means that, when both beams are considered, there will not be a net force of zero 
on a given particle, and so for this imposed uniform spacing the system will not be in 
equilibrium . 



2 

bo 
S 

'o 
a 
ft 

1/5 

0.5 
1 



0.5 

1 2 3 4 5 6 7 8 
Particle number 

Fig. 4. Example of how a force profile like those in Figure [5] (arbitrary units) can be made 
symmetric by altering inter-particle spacings. Now, in contrast to Figure [3] when both 
beams are considered there will be a net force of zero on each particle, and so the sys- 
tem will be in equilibrium. 
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Fig. 5. Collapse of a longer chain, where a shorter chain would be stable. Parameters have 
been selected to give an extreme case where a two-particle chain is stable but a three- 
particle chain is not (parameters as 1161 , but with 1.9 /im diameter spheres). The curves 
"1 of 2" and "2 of 2" show the forces due to the right-going beam for the particles in a 
two-particle chain. Since the force on the second particle is greater than the force on the 
first particle over a range of around 5-12 flm inter-particle spacing, there is a stable trapped 
configuration with an inter-particle spacing of about 12 flm (indicated with an arrow) when 
net effect of both beams is considered. The curves "1 of 3" and "3 of 3" show the same 
forces on the end particles of a three-particle chain. Since there is no spacing for which 
the force on the third particle is greater than the force on the first particle, the chain will 
collapse. Even though there is some enhancement of the force on the third particle over 
what it would be in the absence of the other particles, this is not enough to overcome the 
compressive force due to the effect of the beam on the first particle. 



Finally, Figure [5] shows how a chain above a critical length can collapse. The figure shows 
how the forces on the end particles in a short chain vary with inter-particle spacing for a par- 
ticular set of experimental parameters (different to those used earlier, and carefully selected 
so the collapse occurs at an unusually short chain length). In this case a two-particle chain is 
supported, but if a third particle is added the chain will collapse until the spheres are in contact 
(an effect mentioned in lfl6l ). As pointed out earlier, the compressive force on the first particle 
in the chain is greater if the chain has more particles in it (since that first particle is closer to the 
beam waist). We argued earlier that this would cause the inter-particle spacing to fall until the 
inter-particle repulsive force was increased enough to for the two forces to balance. However, 
for the particular parameters in this figure, at short inter-particle spacings there is in fact an 
attractive force between neighboring particles, and so the three-particle chain collapses once 
the compressive forces have pushed the particles close enough to enter this regime. Our im- 
plicit assumption that the inter-particle light forces are repulsive (it was assumed that radiation 
pressure will dominate) has broken down; near-field gradient force effects have come into play, 
producing a net attractive force between the spheres at close ranges. There is no repulsive force 
to support the chain, and it collapses. 

4. Conclusions 

We have explained the mechanisms behind the formation of optically bound particle chains in 
counter-propagating Gaussian beam traps. The optical binding effect results from the balancing 



of repulsive effects from the light from one particle incident on the next particle in the chain 
and compressive effects due to the background trapping potential formed by the beams (it could 
in fact be argued that this effect is not optical "binding" in the strict sense shown in early 
experiments ifTTl . since the interaction here is largely a repulsive one, with stable chains only 
being formed due to the background harmonic potential of the trap). Here we have used a very 
simple model to successfully explain the trends of closer spacings as more particles are added 
to the chain, and of closer spacings in the center of a chain compared to near its edges. 

While our simple model does not claim to agree precisely with experimental results and with 
the theoretical gold standard of Mie scattering calculations which we have also used (and has 
a number of parameters which must be tuned by hand), there is good qualitative agreement 
between them across a range of model parameters. From this we can conclude that, while 
there is some influence from more sophisticated effects which can only be encapsulated in a 
full vector model based on rigorous solution of Maxwell's equations (such as Mie scattering), 
many of the properties of the trapped particle chains can be understood in terms of a simple 
scalar model. This model can offer strong conceptual insights into the physical mechanisms 
which lead to the observed behavior, which had not previously been fully explained. As well as 
explaining how the inter-particle spacings are regulated, it explains the trend for closer spacings 
with larger N, and the wider spacing close to either end of the chain. It is only through a simple 
model such as the one we have presented that the various complex effects in the experiment 
can be decoupled in order to understand why optical binding occurs under these experimental 
conditions. 
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